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In order to study whether the inter-band nesting can favor superconductivity arising from 
electron-electron repulsion in a three-dimensional system, we have looked at the repulsive Hub- 
bard model on a stack of honeycomb (i.e., non-Bravais) lattices with the FLEX method, partly 
motivated by the superconductivity observed in MgB2. By systematically changing the shape of 
Fermi surface with varied band filling n and the third-direction hopping, we have found that the 
pair scattering across the two-bands is indeed found to give rise to gap functions that change sign 
across the bands and behave as an s- or d-wave within each band. This implies (a) the electron 
repulsion can assist gapful pairing when a phonon-mechanism pairing exists, and (b) the electron 
repulsion alone, when strong enough, can give rise to a d-wave-like pairing, which should be, for a 
group-theoretic reason, a time-reversal broken d-|-id with point nodes in the gap. 



PACS numbers: 74.20.Mn 



I. INTRODUCTION 

Recent discovery of the superconductivity in MgB20 
with relatively high transition temperature (Tc ~ 39K) 
has invoked renewed interests in sp-bonded materials. 
Electronically, the system is a tt electron system on lay- 
ered honeycomb lattice, which immediately reminds us of 
grap^^e intercalation compounds (GIC) such as LiCga or 

While the GIC is considered to be a conventional su- 
perconductor with Tc < 5K, MgB2 has an unusually high 
Tc for sp-bonded materials (with a recent exception of 
Ceo-FET structureQ). Recently, Choi et alB have used 
an ab-initio pseudopotential density functional theory to 
solve Eliashberg's equation numerically, and have renra- 
duced Tc ~ 39K, isotope-effect exponent as ~ O.SQfi, 
and have obtained a gapful BCS pairing, which is cod^- 
tent with experimental results such as-|Specific heatuO, 
tunneJixig and photoemission spectraliinU, penetration 
depthEj, and the Raman spectrata . Thus GIC and MgB2 
both seem to be mainly phonon-mediated superconduc- 
tors. However, to realize a high Tc, electron repulsion 
should not stand in the way of the phonon mechanism, 
so the question may be paraphrased: can the electron 
repulsion stand away from or possibly even assist the 
phonon-mediated pairing. 

On a more positive side, superconductivity from 
electron-electron repulsion itself is fascinating in many 
ways, but there are many open questions. While there 
is a growing consensus that high-Tc cuprates may be 
related to the electron correlation, we are only begin- 
ning to understand the link between the underlying band 
structure and the way in "jsihich the electron-mechanism 
superconductivity appearalJ. Indeed, the way in which 
the superconductivity occurs is sensitively affected by the 
shape of the Fermi surface. Recently, two of the present 



authors proposedtS that multi-band systems should open 
a new possibility of much higher Tc, where a fully gapped 
BCS gap function can appear when the Fermi surface 
consists of disconnected pieces, while the usual wisdom 
dictates that the repulsion-originated superconductivity 
should have, as in the cuprates, a strongly anisotropic gap 
function with nodes. They have conceived and demon- 
strated that, for some two-dimensional lattice models, 
the pair scattering (the Coulombic matrix elements that 
scatter pairs of electrons across the Fermi surface) can 
occur across the pockets, which gives rise to a gap func- 
tion that changes sign across the pockets with the same 
sign within each pocket. 

So the purpose of the present paper is a combination 
of the above two motivations. Namely, we study whether 
the inter-band nesting can favor superconductivity aris- 
ing from electron repulsion in a three-dimensional sys- 
tem, by taking the repulsive Hubbard model on a stack 
of honeycomb lattices as a prototype. There, the honey- 
comb, a typical non-Bravais lattice, provides two pieces 
of the Fermi surface arising from the two bands, while 
the stacking can provide a natural nesting along that di- 
rection. So, if the pair scattering across the two bands 
along the nesting vector works favorably, we can expect 
a pairing from the interband nesting with the BCS gap 
with opposite signs across the two bands. 

So we have studied the systematic dependence of 
the solution of Eliashberg's equation with the multi- 
band fluctuation exchange approximation on the shape 
of Fermi surface by changing the band filling and the 
third-direction hopping. The presence of a strong inter- 
band nesting is found to indeed give rise to gap functions 
that change sign across the two bands, but, if we turn to 
the symmetry within each band, there exist two, nearly 
degenerate modes that behave respectively like s- and 
d-waves. This implies (a) the electron repulsion can as- 
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sist an s-wave like pairing when the phonon-mechanism 
pairing exists as a dominant mechanism. We have fur- 
ther found that (b) when the electron repulsion is strong 
enough, the d-wave is realized, which should be, for 
a group-theoretic reason, a time-reversal broken linear 
combination of two symmetries with point nodes in the 

gap. |-n 

According to the band calculationlj, the Fermi surface 
of MgB2 consists of two tubular networks having a boron- 
2p7r character along with two cylinders of the boro n- Zpg 
characters. As has been stressed by several authorst£lEJ, 
the nesting between these two tt bands are quite good. 
Although pa bands are considered to be important in 
that the electron-phonon interaction is much stronger in 
this band, here we concentrate on the pn bands in or- 
der to focus on the effects of the Fermi surface nesting. 
The GIC, on the other hand, has a much smaller inter- 
layer transfer energy (tz), so that the Fermi surface is a 
cylinder of carbon tt character with no dominant nesting 
vectors. So in the latter part of the paper we shall cover 
both MgB2- and GlC-situations in a systematic variation 
of tz and band filling. The systematic study also serves 
to explore how the inter- and intra-band nesting compete 
in realizing superconductivity. 



The noninteracting band dispersion for 3D honeycomb 



e(k) = 2tz coskz ± 3 + 2 [cosfc^r + cosky + cos{kx + ky)]. 

(2) 

In the two-band FLExElE, Green's function G, self- 
energy S, spin susceptibility x, and the gap function 
(j) all become 2x2 matrices, such as Gap{k), where 
k = (k, iujn) with cj„ = (2?! — l)7rT being the Matsubara 
frequency for fermions. 

The self-energy is given by 



T 



(3) 



where the fluctuation-exchange interaction V^^^ (q) is 



Q/3 



l + C/x'"(g) 



al3 

(4) 



with 



II. FORMULATION 

Now let us start with the case corresponding to 
MgB2. Spin- fluctuation- mediated superconductivity is 
here studied with the fluctuation exchaatte. approxima- 
tion (FLEX) developed by Bickers et alMM. The model 
is a 3D two-band Hubbard model with the repulsion U 
on layered honeycomb lattice, 

A,B 



(1) 



+t^EE<T<i' 



where U is the Hubbard repulsion. The essential ingre- 
dient here is the non-Bravais lattice having A and B 
sublattices, which has two bands within a layer. Here 
we assume a non-staggered layer stacking as realized in 
MgB2 and GIC (Fig.|l|) , which also depicts the inter-layer 
hopping tz- Hereafter we take the intra-layer hopping 
t = -1. 
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FIG. 1. 3D layered honeycomb lattice (the left panel), 
which is topologically equivalent to the lattice in the right 
panel. A and B indicate sublattices. 
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^G„p(k + q)G0^{k). 



(5) 



Here we denote q = (q, ie;) with e; = 27rZT being the 
Matsubara frequency for bosons, and N the number of 
k-points on a mesh. 

With Dyson's equation, 



[G(fc)-i]^ [G"(fc)-i] -^E^^fc), 



(6) 



where G° is the bare Green's function G°g(fc) = [{iujn 



af}\ 

jjL — with the bare energy, we have solved 

eqns.(S)-(p[) self-consistently. 

Tc may be obtained from Eliashberg's equation (for the 
spin-singlet pairing). 



Mapik) 



T 

'n 



V^l\k - k')Gaa'{k')Gi30,{-k')c^a,f,,{k'), (7) 



EE 
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where cj) is the gap function, and the pairing interaction 
V^(^^(fc) is given as 



yrAq)-^u- 



x"'{q) 



l-Ux'"{q) 

USafj. 



X'"{q) 



l + U^"{q) 



(8) 



Tc is determined as the temperature at which maximum 
eigenvalue A becomes unity. 
The susceptibility. 
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Xa/3(k,0) = 



X'"(k,0) 



1 _ [/x"'i (k,0) 



(9) 



Q/3 



may be expressed as diagonalized components, 
XAA + Xbb 



± 



XAA - XBB 



IXABi 



(10) 



Throughout this study, we take N = 32^ k-point meshes, 
and the Matsubara frequencies ojn from —{2Nc — l)7rT 
to {2Nc — l)TrT with Nf. = 4096, which gave converged 
results. 



III. RESULT 
A. MgBa 



which breaks the time-reversal symmetry. This combina- 
tion has point nodes on the Fermi surface, which is curi- 
ous but natural as evident from Fig.^ which superposes 
two sets of nodal planes to show how the nodal planes 
intersect each other along some lines for the doubly de- 
generate function and how these lines in turn intersect 
the closed Fermi surface. 



The reason why 's-wave' is only subdominant may be 
traced back to the Fermi surface, which is rather ex- 
tended in k space in this particular case, so that there 
is an appreciable contribution from the intra-band pair 
scattering that the d-wave can exploit. 



Let us first discuss the spin structure. We have fitted 
the shape of thcj-fiermi surface to that obtained by the 
band calculationt2l to have — 0.65, n — 1.03, U = 1.5 
at T p- 0.01, where n is band filling (n = 1 for half 
filling)La. The Fermi surface, Fig.^, consists of two sets 
of tubular networks corresponding to the bonding and 
anti-bonding n bands. The spin-susceptibility x+(k, 0) 
displayed in the same figure shows a sharp peak around 
(0,0, tt) reflecting a good nesting along the z direction 
(arrows in Fig||). 

Figure ^ shows the gap function obtained from Eliash- 
berg's equation. The solutions having the largest A are 
gapless (/>di and (/>j_ which are degenerate. As expected 
from Kuroki-Aritallfl, a gapful (j)s does exist, although its 
A is slightly smaller than that for the d-wave. We have 
called the former solutions 'd-wave' in that the gap func- 
tion changes sign as -I 1 — azimuthally (i.e., within 

each band), while the latter gap 's-wave' in that it does 
not. 

As a hallmark that the scattering of the (intra-band) 
pairs across the inter-band nesting is exploited, all these 
solutions indeed have 

0AA(fc)0BB(fc') < 0. (11) 

We can confirm in Eliashberg's eq.(^ that the sign 
change across the two bands works favorably in increasing 
A if we note V^^ {k — k') > (peaked around k — k' — 
(0,0, tt)) and a relation for multiband Green's function 
Gc.pi-k) = Glfiik) which gives GAB(fc')Gls(-fc') > 0. 

The fact that the dominant 'd-wave' solutions are dou- 
bly degenerate can be understood by a group theoretical 
argumentnJ, in which these solutions belong to Tg rep- 
resentation for the honeycomb system (while 's-wave' to 
T^). The true gap function below Tc to maximize the 
gap should be a linear combination of the two d-waves, 

(t>di + i'l>d2, (12) 




FIG. 2. Fermi surface (left panel, blue: bonding band, red: 
anti-bonding band) and x+ (right panel) with — n for 
U = 1.5, n = 1.03, = 0.65, T = 0.01. Yellow arrows indi- 
cate the nesting vector. 
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FIG. 



The 



sign 



of the gap functions (jifj^ (top), (jiiti^ (middle), and (jis (bottom) 
for U = 1.5, n = 1.03, = 0.65, T = 0.01. We have displayed 
the sign of the gap on the Fermi surface (left panels: bonding 
band, right panels: anti-bonding band) along with the nodal 
planes displayed in green. 
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FIG. 4. Left panel superposes and (I)d2 for the bonding 



band. Right panel plots l^^i + i4>d2 
circles denote point nodes in cpai + ii 



for fcz 



0. The red 



If we turn to the temperature dependence of A in FigJ|, 
we can see that A is significantly smaller than unity even 
for T 0.01 This implies that the spin fluctuation 
alone is not strong enough to realize the d-wave supercon- 
ductivity in this temperature range for the band filling 
and t, taken here. 



On the other hand, the A for the gapful 's-wave' pair- 
ing is seen to have nearly the same magnitude as that of 
the gapless d-wave, although A is again small. Thanks to 
the absence of nodes on the Fermi surface, this one has a 
gapful pairing, which is eligible for assisting the phonon- 
mediated pairing if the electron-phonon intcractionjjs 
considered on top of the electron-electron interaction^^. 
So we conclude that inter-band spin fluctuations can 
work cooperatively with intra-band phonons to realize 
a gapful superconductivity. 

U=1.5 n=1.03 tz=0.65 




o.oe 



FIG. 5. The eigenvalue of Eliashberg's equation A versus 
temperature. The largest solutions are doubly degenerate. 
Lines in this and following figures are guide to the eye. 



B. Optimization for the 3D honeycomb lattice 



Let us depart from MgB2 to move on to the strong cou- 
pling regime in search of superconductivity from electron 
repulsion alone. In general, pairing instability mediated 
by spin fluctuations in 3fi.rSystems is definitely weaker 
than that in 2D systemsllIiE2l. This has been shown in 
a FLEX study fyr the Hubbard model by three of the 
present authorgl3, who have identified its reason in the 
k space volume fraction of the effectively attractive pair 
scattering region that is much smaller in 3D than in 2D. 
Furthermore, 3D systems have a strong tendency toward 
various magnetic orders, so that to identify the 3D sys- 
tems that favor superconductivity from electron repul- 
sion becomes a challenging problem. 

Here we have optimized the pairing instability in the 
layered honeycomb lattice by varying the intcr-laycr hop- 
ping tz and the band filling n. Namely, we have searched 
for sets of parameter values that give large values of A 
without encountering antiferromagnetic instability at low 
temperatures. The resulting best parameter set is found 
to be [/ = 8, n = 1.15, t^ = 0.7 (inset of Fig.|), for which 
the Fermi surface and the (inverse) spin susceptibility are 
shown in Fig.^. In Fig.^ we can see that the estimated 
Tc ~ 0.001 (Aniax ~* !)■ The pairing symmetry is again 
0di + i4>d2 ■ 
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FIG. 6. Fermi surface (left) and x+(^z = (right) for 
the optimized parameter set of f/ — 8,n — 1.15, = 0.7 at 
T = 0.01, where the pairing symmetry is d-wave. 



U=8.0 n=1.15 tz=0.7 
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FIG. 7. A versus temperature for the optimized parameter 
set employed in the previous figure. Dotted line is a spline 
extrapolation to lower temperatures. 



C. Effect of the inter-band nesting 

Apart from the above optimization, we have system- 
atically explored how the inter-band nesting affects the 
pairing in the 3D honeycomb lattice. Since we wish to 
separate out the effect of the band filling and the hop- 
ping in the z direction, we have done this along two paths 
displayed in FigJ|. The first path starts from the param- 
eter values corresponding to the p7r bands in MgB2 , while 
the second path includes the parameter regime {t^ < 0.2) 
which corresponds to GIC except for the value of U. The 
result along the first path is shown in Fig.^, where the 
n-dependence of Amax is plotted for U = 1.5, tz = 0.65 
at r = 0.01. We can see that the pairing instability be- 
comes weaker as n is increased, which is natural since the 
inter-band nesting becomes degraded along this path. 

The result along the second path is displayed in Fig.[ic|, 
where A is plotted as a function of for n = 1.2 at 
T = 0.01. Here we have adopted a rather large U = 8, 
because we want to have sizeable A over a wide range of 
tz including the case of bad nesting. There, we have cov- 
ered both the layered case {t^ < 1) and a quasi-lD case 
(tz > 1). As indicated in Fig.^ antiferromagnetism 



occurs (i.e., x ~* oo) when the inter-band nesting be- 
comes too strong as — > 1 both from below and from 
above. Before the transition to antiferromagnetism oc- 
curs A increases both from below and from above, indi- 
cating that inter-band 3D nesting is effective, the layered 
case {tz — 0.8) is more advantageous than the case of 2D 
{tz = 0) and quasi-lD {tz > 1) on the present path. 



tz\ 


1.03 1.05 1.1 1.15 1.2 1.25 


2.3 




U=8,0 


0.65 


U=1.5 






► 










FIG. 8. The paths we have focused on in the present study 
in the parameter space of the inter-layer hopping and the 
band filling n. 
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FIG. 9. A versus n for [/ = 1.5, = 0.65, T = 0.01. The 
d-wave solutions having the largest A are doubly degenerate. 
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FIG. 10. A versus ior U = 8,n = 1.2, T = 0.01. Grey 
region indicates the antiferromagnetic phase. The insets de- 
pict the shape of Fermi surface at five points on the horizontal 



IV. CONCLUSION 



In conclusion, we have studied the possibility of spin- 
fluctuation mediated superconductivity in 3D honey- 
comb lattice systematically. We have shown that if we 
take the parameter set corresponding to the pir bands 
in MgB2, the spin fluctuation favors the gapful pairing, 
which suggests that the electron correlation can help the 
phonon in forming the Cooper pairing. Experimentally, 
the electron repulsion acting constructively may be con- 
firmed if some phase-sensitive method can detect the gap 
function having opposite signs in two vr-bands. When 
strong enough, the electron repulsion alone will give rise 
to a d-wave pairing, with a time-reversal broken +i(t>d2 
symmetry associated with degenerate representations in 
the non-Bravais lattice with peculiar point nodes on the 
Fermi surface. 
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